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PERFORMANCE OPTIMIZATION OF AN MHO GENERATOR WITH 
PHYSICAL CONSTRAINTS 

C. C. P, Plan, G. R. Seikel and J. Marlin Smith 
NASA Lewis Research Center 
Cleveland, Ohio 4413C 


ABSTRACT 

A technique has been described which optimizes 
the power out of a Faraday MHO generator operating 
under a prescribed set of electrical and magnetic 
constraints. The method does not rely on complicated 
numerical optimization techniques. Instead the 
magnetic field and the electrical loading are 
adjusted at each streamwise location such that the 
resultant generator design operates at' the most 
limiting of the cited stress levels. The simplicity 
of the procedure makes it ideal for optimizing 
generator designs for system analysis studies of 
power plants. The resultant locally optimum channel 
designs are, however, not necessarily the global 
optimum designs. 

The results of generator performance 
calculations are presented for an approximately 2000 
MWg size plant. The differences between the 
maximum power generator design and the optimal design 
which maximizes net MHO power are described. The 
sensitivity of the generator performance to the 
various operational parameters are also presented. 

INTRODUCTION 

MHD power generation, using fossil fuel, 
promises both high efficiency and low cost of 
electricity. One of the most critical components in 
the MHD/steam power plant is the MHD generator. Its 
performance strongly influences the overall 
performance of the power plant. The generator design 
should be selected to maximize the overall plant 
efficiency. The optimal design is, however, a 
function of the generator operating parameters 
including the combustor oxidant preheat temperature, 
pressure, and Oxygen enrichment; the magnetic field 
strength and profile; the generator electrical load 
parameter profile and the Mach number; and the state 
of generator technology in terms of allowable 
electrical stresses, etc, 

This paper describes a technique which predicts 
the variations of the magnetic field and the 
electrical loading along a Faraday MHD generator for 
optimum power out. Illustrative results of generator 
performance calculations using this technique are 
presented for an ECAS-2 size plant - 2000 MWg. The 
sensitivity of generator performance to the various 
operational parameters is investigated, 

OPTIMIZATION PROCEDURE 

A prescribed set of electrical and magnetic 
constraints is imposed in the present procedure to 
optimize the power of a generator with specified 
length, diffuser exit pressure, and the streamwise 
distribution of either velocity or Mach number. The 


magnetic field, B, and the electrical loading, K, are 
adjusted in such a way that the generator attempts to 
operate over a large fraction of its length at 
specified stress limits; thus, the generator 
performance is maximized. The electrical "loading is 
defined as: K = jEy, where Ey is the Faraday electric 

uB 

field and u is the channel axial velocity. The 
prescribed constraints may include limits on the 
Faraday current density, Jy, the Hall field. Ex, the 

total electric field, E ^x ^ ^y ’ 

parameter, B . The channel design is such that at 
each streamwise location all of the following 
conditions are satisified: 
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Appearing in equation 12) is the electrical 
conductivity, a. The magnitudes of B and K are 
determined at each station by only two appropriate 
criteria. Which two criteria prevail at a location 
depend on the generator operating conditions and the 
assumed values for the electrical and magnetic 
constraints. The values of B^jax, 

Jycrit> Bcrit should be selected to reflect 

the current or projected states of magnet and channel 
technology. The channel operatring life would also 
be strongly influenced by the values assigned to 
these limits. 

The insulating sidewalls of the Faraday 
generators are subjected to both axial and transverse 
electric fields. In order to limit this electrical 
stress on the sidewalls, the total electric field 
strength is bounded by condition (4) in the present 
procedure. When selecting a suitable value for 
^criti oils ™^st consider such factors as: the 

particular sidewall design, the slag layer behavior, 
and the wall temperature. Whether the transverse 
electric field, Ey, should be limited instead of E 
would also depend on these factors. In that case, 

Ey < Evcrit replace condition (4). 

Excessively high values of Ex should be 
avoided within the generator in order to prevent 
interelectrode breakdowns. This critical Ex has 
been verified experimentally to be approximately 
4 kV/m. The traditional approach to generator design 
is to limit Ex by limiting of the generator 
magnetic field via condition (1). In the present 
formulation, one has the option of using either 
condition (3) or condition II) to accomplish this. 

The bulk plasma non-uniformities in 
combustion-driven generators can become severe at 
high values of Hall parameter as in inert-gas 
generators. This effect, however, genera’ly occurs 
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in open-cycle generators at a much higher value of 
3. One can conveniently use condition (1) to limit 
the Hall effect. However, electrothermal and/or 
magnetoacoustic wave instability in combustion 
generators are still poorly understood and the 
critical 8 associated with these wave phenomena have 
yet to be defined; a Scn't value of 4 is widely 
assumed as a reasonably safe value. 

Knn'p is the minimum value of the load factor 
and is used as a control parameter in the calculation 
procedure. For a given generator inlet condition, 
Kmin is adjusted in an iterative manner until the 
desired channel length and exit pressure are 
satisfied. The calculated MHD power for the various 
operating conditions are then compared to ascertain 
the optimum MHD power condition. 

An optimum magnetic field profile is calculated 
by the present procedure. This can then be used as 
guidelines by the magnet designer. After the magnet 
design is finalized, the generator loading should be 
reoptimized using the actual magnetic field profile 
in place of conditions (!) and (6). 

An approach to the optimization of the Faraday 
channel loading to maximize the gross MHD power has 
been previously suggested by Doss and Geyerl. 

Their sectional optimization technique uses the Hall 
field strength and the Faraday current density as 
criteria for controlling the load factor 
distribution. The magnetic field was kept constant 
at a value of 5 Tesla. The controlling criteria to 
maximize the MHD power in the present formulation 
become identical to those of reference 1 if 
B = 5 Tesla replaces condition (6), ^crif^ “> 

Ecrit However, as discussed below, the 

maximum MHD power design is no optimum from a system 
viewpoint. 

CHANNEL PERFORMANCE CALCULATIONS 

To illustrate results of generator performance 
calculations using the technique described, generator 
design was examined for an ECAS-2 size plant^. The 
fuel is Illinois #6 coal and the seed is potassium 
carbonate. Other design parameters are listed in 
Table I. Boundary conditions are specified at the 
two ends of the MHD generator. The total enthalpy at 
the generator inlet is given as the enthalpy after 
the combustion less the heat loss in the combustor 
and nozzle. The ECAS-2 diffuser exit stagnation 
pressure of 1.14 atm, diffuser pressure recovery 
coefficient of 0.7, and generator length of 25 meters 
are assumed in this study. The generators are also 
assumed to be operating at approximately constant 
Mach number (■YgM^ = constant). The comparison 
between the performance of the approximately constant 
Mach number and the constant velocity designs and the 
sensitivity of generator performance to shorter 
channels and less efficient diffusers are presented 
elsewhere2>3. 

The generator calculations are performed with a 
quasi -one-dimensional flow model. The model consists 
of an inviscid central core flow with boundary layer 
developing along the walls.;. The turbulent boundary 
layers are treated with a momentum integral method. 
The combustion chamber conditions and the 
thermodynamic and transport properties of the 
combustion gas are calculated following Svehla and 
McBride'^. The electrical properties of the gas 
were computed as in reference 5. The electrode 


voltage drop distributions are taken to be of the 
form: V(j = a + b 6* where 6* is the displacement 

thickness. The values of a and b should generally be 
functions of the generator operating conditions; 
however, they are kept constant in the analysis to 
limit the complexity of the generator model. The 
value of a = 100 V and b = 20.14 V/cm are taken from 
reference 6. 

The t^ical variation of power out for an MHD 
generator is shown in figure 1. The gross MHD power 
as a function of the minimum load factor, is 

shown in figure la. The net channel power out, MHD 
minus compressor power, versus the MHD combustor 
pressure. Pc, is shown in figure lb. In the 
computation of compressor power, an initial oxidant 
temperature of 60F and a compressor polytropic 
efficiency of 0.898 are assumed. A pressure drop of 
0.163 Pq is assumed between the air compressors and 
the MHD combustor. As indicated in figure la, there 
is a maximum power and corresponding combustor 
pressure which the specified length Mn this case 25 
meters) generator can be operated at without 
exceeding the specified electric and magnetic 
stresses. Or in effect, for a given combustor 
pressure there is a minimum length channel. In the 
majority of the cases, this minimum length condition 
coincides with the maximum power out condition for 
the given chamber pressure. The optimum operating 
combustor pressure is, however, at the maximum of the 
net power curve of figure lb. 

The channel design is vastly different for the 
maximum power and maximum net power generators. In 
figure 2, distributions of the electrical variables 
and the magnetic field for these two generator 
designs are compared. For this example, the maximum 
net power generator operates at a combustor pressure 
of 9 atm. The channel has a constant load parameter, 
K= 0.798. Conditions 14) and (5) dominated over a 
major portion of this generator. The maximum power 
generator has a higher combustor pressure of 12.786 
atm and variable loading is required. The load 
factor tends to decrease near the channel entrance 
(trend dictated by conditions (2) and (4)), then 
increases rapidly downstream (trend dictated by 
conditions (3) and (4)). The Bcrit condition 
prevails beyond - 23 meters causing the drop in the 
magnetic field and K profiles. The Jy, E^^, and K 
profiles for the maximum power design'^are similar to 
the results of Doss and Geyerl. Jhe optimum net 
power generator design is preferred, however, from 
overall plant efficiency considerations2,3. 

Typical net generator power variations are 
plotted in figure 3 for different oxidant preheat 
temperatures. The trend towards higher plant 
operating pressures with increasing' preheat 
temperatures is clearly shown. The thermal input to 
the generator was maintained at the 5373 MWth level 
in these calculations by adjusting the mass flow rate 
for the various preheat temperatures. 

The sensitivity of generator net power to the 
design constraints is shown in figure 4. Figure 4a 
shows the net power versus chamber pressure for 
different sets of assum>.d electric and magnetic field 
limits. For E^rit " ^ ^Xcrit ” 2.5 kV/m, 

and = 1 A/cm2 (approximately the 

electrical stress limits in present generator 
endurance tests), no further increase in power out is 
possible if the peak magnetic field strength is 
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increased (except at nonoptimally high values of 
combustor pressure). Only by increasing the limits 
on the electrical constraints can the higher magnetic 
field be utilized. As shown, however, operation at 
higher values of electrical stress offers significant 
performance improvement. 

Figure 4b shows the sensitivity of net generator 
power to the assumed value of The 

minimal effect of od net power is due to 

the fact that very small part (if any) of the maximum 
net power generator operates at the 

limit. By contrast, a large portion of the maximum 
MHD power generator is stressed to the allowable Jy 
limit. For the sample calculations of figure 4b, the 
net power decreases by only - 1.5;^ when the value of 

“^ycrit lowered to .5 A/cm^. Operation at 

the reduced current densities might greatly benefit 
the channel reliability and lifetime. 

The influence of Mach number on channel 
performance is shown in figure 5. Results are 
presented for different and B^ax 

limits. Variable pressure recovery coefficient as a 
function of Mach number was assumed for the 
supersonic diffusers. The optimum Mach number 
decreases as the peak magnetic field is increased. A 
slightly higher Mach number is desirable for 
generators operating at the higher current 
densities. The reason for the insensitiveness of net 
power to is as mentioned previously. 

Enthalpy extraction for supersonic channels are 
substantially lower than for the subsonic channels 
operating at the same electrical stress levels. This 
results from the need in the supersonic generators to 
restrict the magnetic field strength in order to 
satisfy the given constraints. 

CONCLUDING REMARKS 

The "local" optimization approach of the present 
formulation or that of reference 1 may not yield a 
global optimum design, as discussed in reference 7. 
Whether a global optimization with constraints on 
Ey or E will still give better maximum net P^hd 
than the present technique needs to be examined. 

Under the present formulation, the maximum net power 
channel design is forced to a constant generator 
loading by condition (5). Whether this is a 
desirable trait needs to be verified by comparing 
with a global optimization. 
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Table I - Major Parameters 


Coal type 111 inois #6 

Moisture content of coal delivered 2 

to combustor, percent 

Oxidizer preheat temperature, F 2000, 2500, 

3000 

Combustor pressure, atm variable 

Combustor fuel-oxidizer ratio 1.07 

relative to stoichiometric 

Combustor slag rejection, percent 85 

Generator type Faraday 

Potassium seed, seed-coal weight .108 

ratio 

Diffuser exit pressure, atm 1.14 

Diffuser pressure recovery 0.7 

coefficient 

Generator length, meters 25 

Compressor poTytropic efficiency 0.898 

Thermal input to combustor, MW^-j, 5373 
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COMBUSTOR PRESSURE, /-LOAD LINE FOR 
atm / 25 m CHANNELS 



*^min 


(a) GENERATOR POWER AS A FUNCTION OF MINIMUIVl 
LOAD PARAMETER, 



(b) MHD GENERATOR POWER MINUS COMPRESSOR 
POWER AS A FUNCTION OF COMBUSTOR PRES- 
SURE. P^. ^ 


Figure 1. - Typical MHD generator power variation. 
Oxidizer 2500° F preheated air, generator length 
25 m, thermal input 5373 MW, ygM^ = const. , 

*'^inlet “ ^crit 

V/m, =• 10 000 A/m^, = 4. 


r«"' 


MAX. NET POWER 

MAX. POWER 


(a) GENERATOR LOAD PARAMETER VARIATION. 


(b GENERATOR MAGNETIC FIELD VARIATION. 


STREAMWISE DISTANCE, m 
(c) ELECTRODE CURRENT DENSITY VARIATION. 

Figure 2. •- Parameter distributions for typical MHD gen- 
erator ease. Oxidizer 2500° F preheated air, generator 
length 25 m, thermal input 5373 MW, y^M^ = const., 
Mjniet = 0- 9, Ecnt = 4000 V/m, E„ = 2500 V/m, 

J . = 10 000 A/m2, Pcrit = 4. 
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COMBUSTOR PRESSURE, atm 

Figure 3. - Net MHD power as a function of combustor 
pressure for various oxidant preheat temperatures. 
Generator length 25 m, thermal input 5373 MW, 
yjM^ = const. , Mj^iet = Ecrit " 

E„ . = 2500 V/m, Jw . = 10 000 A/m^, = 4. 
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(a) SENSITIVITY TO ELECTRICAL STRESS AND MAGNET 
ASSUMPTIONS. Jy = 1 A/cm^, Pcrit = 4. 
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Figure 5. - Maximum net power as a function of gener- 
ator Mach number. Oxidizer 2500° F preheated air, 
generator length 25 m, thermal input 5373 MW, 
const. , E(,pj{ = 4000 V/m, = 2500 V/m, Pj^rit " 

diffuser pressure recovery coefficient = 0.7 for 
Mj^l; =0.98 -0.28 Mj for Mj^l. 



